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ABSTRACT 

In this paper we use a path-integral approach to represent the Lyapunov exponents of 
both deterministic and stochastic dynamical systems. In both cases the relevant correlation 
functions are obtained from a (one-dimensional) supersymmetric field theory whose Hamil- 
tonian, in the deterministic case, coincides with the Lie-derivative of the associated Hamil- 
tonian flow. The generalized Lyapunov exponents turn out to be related to the partition 
functions of the respective super-Hamiltonian restricted to the spaces of fixed form-degree. 



Published in: Chaos, Solitons and Facials, Vol.4 no 7 (1994) 1117 



1. INTRODUCTION 



During the past two decades a lot of efforts have been devoted to the study of the so-called 
"stochastic" properties of deterministic dynamical systems. Most of the work concentrated 
on finding appropriate order parameters which could be used to classify dynamical systems 
according to their degree of stochasticity' 11 P1 . Rigorous definitions of integrable, KAM-, 
ergodic, weakly mixing, mixing, C-systems, etc. were given and it was established that 
this ordering of the systems amounts to an increasing "chaoticity" or "stochasticity" of the 
motion. It also became clear that some of the stochasticity properties of the system (like 
ergodicity) were encoded in the spectrum of its Liouville operator. 

One of the order parameters which has been studied most extensively is the Kolmogoroff- 

[3] 

Sinai (KS) entropy and similar entropy-like quantities . A priori, the KS entropy is defined 
in a rather abstract manner*. It has been shown,' 41 however, that it is related to the sum 
of the positive Lyapunov exponents W . These exponents, loosly speaking, are a measure of 
the rate at which nearby trajectories fly away from each other. For deterministic systems 
the Lyapunov exponents are computed from the properties of a single trajectory, i.e., they 
are labelled by the initial point of the respective trajectory. More recently the concept of 
Lyapunov exponents has been generalized also to stochastic systems' 61 . In this case the 
Lyapunov exponents, like all observables, are obtained by averaging over (infinitely) many 
trajectories. This averaging can be done in different, inequivalent ways [see ref.9]. It is well 

ryi 

known that stochastic systems with a Langevin dynamics can be formulated via path- 
integrals and, once the noise is integrated away, they are equivalent to a one-dimensional 
super symmetric field theory. As was pointed out before' 8 " 91 there exists a special class of 
generalized Lyapunov exponents which acquires a natural interpretation in terms of this 
super symmetric field theory. In fact they are simply given by the lowest eigenvalues of the 
super-Hamiltonian at fixed fermion number. Hence they can be deduced from the asymptotic 
behaviour of the corresponding partition function. 

The relation between Lyapunov exponents and supersymmetry is much more general 
than it might appear from the previous work on stochastic systems. In fact, let us consider 
a trajectory <p(t) of some dynamical system, and let us visualize <f>(t) as a "bosonic field" 
in a one-dimensional field theory with t and <fi parametrizing the "spacetime" and the 
"target-space", respectively. Let us assume that the lagrangian £ of the system is invariant 

* See ref.[3] for the definition and the original references. 
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under a supersymmetry transformation of the form 



6<j>(t) = ec{t) 



(1.1) 



where c(t) is the anticommuting "super-partner" of (f>(t). Furthermore, assume that C 
leads to an equation of motion for (f) which is of the general form 



where V is some vector field on the "target- space" and rj(t) is any function indepen- 
dent of 0, and which does not change under the supersymmetry transformation. Then 
supersymmetry implies that the "fermionic partner" c(t) evolves according to 



where V is the Jacobi matrix (see below) of V. Obviously the equation of motion of the 
super-partner c(t) is the linearization of the "bosonic" equation of motion (1.2). Hence the 
dynamics of c(t) contains information about the stability properties of the 4>— trajectories. 
In particular, if two nearby trajectories fly away from each other exponentially fast, this will 
manifest itself by exponentially growing eigenmodes of c(t). 

The above argument about the connection between supersymmetry and (generalized) 
Lyapunov exponents is very general. It applies to stochastic and deterministic systems 
alike. However, only for stochastic systems the field theory approach we mentioned has 
been widely used in the literature 18 ' 91 , whereas for deterministic systems the relevant field 
theory was introduced only recently! 101 In the stocastic case the pertinent field theory is the 
euclidean-time version of supersymmetric quantum mechanics (see ref. [16] for details). In 
this context the equation of motion (1.2) is the Langevin equation with a white noise r)(t). 
From the point of view of supersymmetric field theory, eq. (1.2) is the Nicolai mapping (see 
the third of ref. [16]) relating (j) to the Gaussian field rj. In the deterministic case rj(t) = 0, so 
that (1.2) becomes a generic first order evolution equation. We are particularly interested 
in Hamiltonian systems, in which case the vector field V is the symplectic gradient of some 
Hamiltonian H. In ref. [10] we have set up a field theoretic formalism for such systems, in 
particular we introduced a path-integral formulation for deterministic systems. The relevant 
Lagrangian is indeed invariant under a supersymmetry of the form (1.1). 



(1.2) 



c(t) = v'WMt) 



(1.3) 
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The paper is organized as follows. In section 2 we briefly review some basic facts about 
Lyapunov exponents. Then, in section 3, we introduce the path-integral for classical hamil- 
tonian systems. In section 4 we present the observables of this theory which measure the 
(ordinary) Lyapunov exponents refering to a single trajectory. Next we show , in section 5, 
that the partition function of the classical super-hamiltonian is related to a set of general- 
ized Lyapunov exponents which constitute a classical analog of the one found by Graham 
in the stochastic case. In appendix A we give the corresponding discussion for stochastic 
systems. There we generalize previous work for one dimensional configuration space 18 ' 91 to 
higher dimensions in order to elucidate the geometrical structure underlying the higher di- 
mensional Lyapunov exponents. In a second appendix (B) instead we briefly indicate the 
relation between our path-integral and the thermodynamic formalism of Ruelle 1231 with the 
hope to come back in the future to a more complete study of the relations between the two. 



2. BRIEF REVIEW ABOUT LYAPUNOV EXPONENTS 

Let us consider the differential equation 

j t r(t) = h a m)) (2.i) 

where h is a vector field on some N-dimensional manifold M.n with local coordinates <p a . 
(Later on h will become the hamiltonian vector field and 

A4n — M-2n a symplectic manifold.) Let 0o) (the subscript "] c j" is for classical) 

be the solution of (2.1) with initial condition = 0; 0o) — <f>o- We can then define a 

matrix associated to 

s?(t; to) = to) (2-2) 

This matrix is usually known as " Jacobi matrix" [101 and it describes how small changes of 
the initial point affect the solution at time t. The Jacobi matrix is a solution of the linear 
equation 

[d t 5 a b - d b h a ($ cl (t; <f> ))]S b c (t; O ) = (2.3) 

with 5^(0; 0o) = 5%. Eq. (2.3) is also called the "equation of the first variations", in fact if 
we displace the initial point <po by an infinitesimal amount 5<p a (0) then the dispacement 
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at later times is given , to first order, by the Jacobi field 



5<p a (t) = S%(t; 00)5/(0) (2.4) 

Let us pick some smooth Riemannian metric on AAn and let e a denote a unit vector in 
the tangent space T^M. Then A(t; </>, e) = \\S(t; <f>)e\\ is called the coefficient of expansion 
in the direction of e . If 

limsup -In \(t; (/), e) > (2.5) 

there is then an exponential divergence in the direction of e. In ref.[5] ( and summarized in 
the second of ref. [3]) it was shown, under rather general conditions, that: 



i) the one dimensional Lyapunov exponent for the direction of e, as defined by 

A«(e; 0) = limsup -ln\\S(t; <p)e\\ (2.6) 

t^oo t 



exists for every vector e G T^M. 

N 



ii) there exists a basis (called "normal basis" W ) {ej, 1 < i < N} for T^,.M 



such that 

N N 

= in fJ2 xil) ( 2 - 7 ) 



i=i 



where the "infimum" is taken over all possible bases {e^} of T^M.^. It is then easy 
to realize that for any e G T^M, we have 

A«(e;0)G{A«( ei ,0);l<z<7V} 

The numbers A^fe; </>), 1 < i < N, are called the one-dimensional Lyapunov charac- 
teristic numbers or exponents at (j). We shall write 
Aj- 1 ^) = X^\ei, cp) and choose the labelling such that 

aS 1) (0)>a( 1) (0)>..->A«(0) (2.8) 

The numbers {A^|l < i < N} are not necessarily distinct. 
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The higher dimensional Lyapunov exponents X^ are defined in a manner similar to the 
1-dimensional case. Let {e\, • • • , e p } be a set of p < N orthonormal vectors in T^M. 
They span a p-dimensional parallelotope V p = e\ A e2 A • • ■ A e p . The exponents are a 
measure for the exponential growth of the volume of V p : 

X {p \V p ; 0) = limsup -ln\\S(t; <j>)ei A S(t; <p)e 2 A ■ • • A (p)e p \\ (2.9) 

It was shown M that, for almost all initial V p s, A^ = A^ where A^ is given by the sum 
of the p largest one-dimensional Lyapunov exponents 

A^M^Aj 1 ^) (2.10) 

Besides the numbers {X^\<f))}, which are uniquely associated to the flow h a and the point 
4>, there is a quantity that is associated to the flow itself and not to any point in particular. 
It is the Kolmogoroff-Sinai (KS) entropy which we mentioned at the begining and which 
has played a leading role in detecting the transition from ordered to stochastic motion . 
According to a theorem by Piesin W , the KS entropy can be related to the positive Lyapunov 
exponents: 

KS= J #[ Yl A *^)] 
Mn M<P)>o 

This is the main reason why the central objects to study are the Lyaupunov exponents. 

So far we have discussed Lyapunov exponents for deterministic systems only. For stochas- 
tic systems various inequivalent versions of Lyapunov-like quantitites have been studied in 
the literature 161 . Here we only mention the one employed by Benzi et al.' 81 and by Graham 191 . 
For the system (1.1) with 77 a white noise, say, it is not possible to define Lyapunov expo- 
nents for individual trajectories, but only for averages. In ref.[8,9] a generalized Lyapunov 
exponent was defined by 

Al = limsup -In (trS(t)) (2.11) 

where (•) denotes the stochastic average over closed trajectories of length t. Here S(t) is 
the Jacobi matrix evaluated at the end-point of the trajectory (actually in ref.[8,9] only 
systems with one degree of freedom where considered so that S(t) did not have indices). 
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3. THE PATH-INTEGRAL FORMULATION 
OF CLASSICAL HAMILTONIAN DYNAMICS 

In classical mechanics (CM) the propagator P (02,^2 |<^i,^i), which gives the classical 
probability for a particle to be at the point 4>2 at time t 2 , given that it was at the point (pi at 
time t\, is just a delta function 

P(<hM<h, h) =5 2n (<f>2-<s> c i(t 2 , fa)) (3.1) 

where $ c /(t, 4>q) is a solution of Hamilton's equation 

4> a {t) = u ab d b H{(t)(t)), with u ab u bc = 8 a c (3.2) 

subject to the initial conditions 4> a {t\) = 4>\ Here H is the conventional hamiltonian 
of a dynamical system defined on some phase-space M.2n with local coordinates (p a ,a = 
1 • • • In and a constant symplectic structure u> = ^u ab d(f) a A d(j) b . 

The delta function in (3.1) can be rewritten as 

N-l 

8 2n (<f> 2 - <Mt 2 ,0i)) = { II / d ^i) 52n ihi) - ^(^,0o))}^ 2n (02 - <MWi)) (3-3) 

i=i J 

where we have sliced the interval [0,t] in N intervals and labelled the various instants as t{ 
and the fields at U as <p^y Each delta function contained in the product on the RHS of 
(3.3) can be written as: 

In 

^ 2n (<%) - MU, <M) = II S $ a - ^ ab d b H) lt det[5Zd t - ^(uacMdcHte))] ^ (3.4) 

a=l 

where the argument of the determinant is obtained from the functional derivative of the equa- 
tion of motion (3.2) with respect to (j>u\. Introducing anticommuting variables c a and c a to 
exponentiate the determinant, and the commuting auxiliary variables A a to exponentiate 
the delta functions, we can re-write the propagator as a path-integral. 

P{hM\<t>iM) = J VcpVXVcVcexp iS (3.5) 

where S = dt C with 

I = X a [4> a - u ah d b H(4>)\ + ic a {5 a h d t - d b [u ac d c H(<p)})c b (3.6) 

In the path-integral (3.5) we have used the using the slicing (3.3) and then taken the limit 
of N — > oo. This limit has to be taken with some care and some normalization factors might 
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appear in eq. [3.5], but they are of no importance for our discussion. Holding and c both 
fixed at the endpoints of the path-integral, we obtain the kernel' 101 , K((j>2, C2, *2 1 01 , ci, ^l), 
which propagates distributions in the space (0, c) 

£?(02,C2,t 2 ) = y C? 2n Cl K(02,C 2 ,t2|01,Cl,tl)^(0l,Ci,ti) (3.7) 

The distributions £>(</>, c) are finite sums of monomials of the type 

mc) = ^l. a Mc ai --- cap (3-8) 
The kernel K(-|-) is represented by the path-integral 

ti 

K(02,C2,t 2 |0i,ci,ti) = J V<p a V\ a Vc a Vc a expi J dtC (3.9) 

with the boundary conditions a (ti,2) — 4>i2 an< ^ ^(^1,2) — c i2- The function g of 
eq. (3.8) is the classical analogue of a wave function in the Schroedinger picture. It is also 
easy from here to build a classical generating functional Z c \ from which all correlation- 
functions can be derived. It is given by 

Z d = J V(p a (t) VX a (t) Vc a {t) Vc a exp i J dt{£ + source terms} (3.10) 

where the lagrangian can be written as 

C = Xa^ + lC a t b - H (3.11) 

with the "superhamiltonian" given by 

H = X a h a + ic a d b h a c b (3.12) 

and where h a is the hamiltonian vector field 

h a (<p) = u ah d h H(<p) (3.13) 

From the path-integral (3.10) and (3.11) we can see' 101 that the variables (0, A) 
and (c, c) form conjugate pairs satisfying the (^-graded) commutation relations 

[c a ,c„] =St 

The commutators above are defined in precise terms in ref.[10]. Because of these commu- 
tators, in the " Schroedinger-like" picture 1101 , the variables A a and c a are represented by 
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A a = = —id a and by c a = In this way the hamiltonian (3.12) becomes 

H = -il h (3.15) 

where 

lh = h a d a + c\d b h a )-^ (3.16) 

is the Lie derivative operator. Its bosonic part coincides with the Liouvillian L = h a d a , 
which gives the evolution of standard distributions (0) in phase-space: 

dJ 0) (<P,t) = -l h6 {Q \M = Lg(°\cf>,t) (3.17) 

We see from here that our path-integral is nothing else than the path-integral counterpart of 
the operator approach to CM pioneered by Koopman and von Neumann [nI . The next step 
is to understand the meaning of the ghosts c a . From the lagrangian C of (3.11) we see that 
they obey the eq. 

c a (t) = d h h a (<p(t))c h (t) (3.18) 

which is the same equation as the one for the first variations 5<p a derivable from eq. (3.2). 
So we can say that there is a one to one correspondence between ghosts c a {t) and Jacobi 
fields 5(j) a (t). 

c a (t) ~ 5(f) a (t) (3.19) 

Of course we know that the Jacobi field depends on the trajectory <f>(t) we are varying , 
but the same is with the ghosts c a because in solving (3.18) we had to specify the classical 
trajectory <p{t) to insert in h a . In a Schroedinger-like picture, in which the ghosts c a {t) do 
not depend on t, but are only arguments of the g(<f>, c, t), the correspondence of the ghosts is 
not with the Jacobi fields but with the basis of the cotangent space T^M.2n that is usually 
indicated as d(p a . For the details see ref.[10]. With this interpretation of the ghosts, it is then 
easy [u>1 to re-interpret all the Cartan-Calculus on symplectic manifolds and also to identify 
7i with the Lie-derivative of the hamiltonian flow /^appearing in (3.16). Because the ghosts 
c a form a basis of the cotangent space T^M.2n , f?(0, c) may be considered a p-form valued 
field on M.2n- The Lie-derivative acts then on the components of g (3.8) in the standard 
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manner 



k&,..a P = h%£.. ap + £ d aj h b $l ajibaj+i „. ap (3.20) 

i=i 

The kernel K obeys the Schroedinger-like equation idtK = HK. Consequently the 
time evolution of g is governed by the equation 

idto(<i>,c,t)=HQ(<l>,c,t) (3.21) 

or 

dtoil.a p {M = -ho£la p (M (3-22) 
The explicit evaluation 1101 [12] of the path-integral (3.9) yields 

K(^ > c a> t 2 |0i J ci J ti)=^ 2n )(0§-$S(t 2 ;0i))(5( 2n )(^-CS(fe;ci^i)) (3.23) 
Here and C a d (t) are solutions of the classical equations of motion resulting from C, 

r(t) = h a (0(t)) = u ab d b H{(f>{t)) (3.24) 

c a (t) = <9 & /i a ((/)(t))c 6 (t) (3.25) 

with the boundary conditions $^(ti;0i) = 0f and c\, </>i) = cf. 

The third argument of C5(i; c\, </>i) indicates that is the Jacobi field for the clas- 

sical trajectory $^(t) emanating from the initial point <f>\. Eq. (3.25) is solved by 

c a (t) = SS(t)<*(ti) (3.26) 

if the Jacobi matrix obeys the differential equation 

[d t S a b -MZ(t)]S b c (t) = (3.27) 

with the initial condition S^(ti) = 5% and where 

M£{t) ee d h h a m)) = u ac d b d c H(<P(t)) (3.28) 

The formal solution to eq. (3.27) reads 

t 

S(t) = fexp J dt' M(t') (3.29) 

where T denotes the time-ordering operator. Note that S(t) functionally depends on the 
trajectory <p(t) chosen. Since the latter is uniquely characterized by its initial point <pi we 
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shall write S(t) = S(t; (pi) for the Jacobi matrix of the trajectory emanating from <p\. The 
function S(t;<pi) defines what is called' 131 a "multiplicative cocycle": 

S£(t + r; O ) = S a c (t; $ d (r; 0o))5 6 c (r; O ) (3.30) 

It is well known' 14 ' that S is symplectic*, 5 G Sp(2n), i.e., 

Sc^abSd = u cd 

As a consequence we have that detS = 1 . Using (3.26) we write for the Kernel (3.23) 

K(0 2 , c 2 , t 2 \4>i, ci, ti) = 5^\<p a 2 ~ Kl(t2; <Pi))^ 2n \c a 2 - SS(t 2 ; <Pi)c\) (3.31) 
In what follows we shall frequently exploit the fact that K is normalized, 

J rf 2n 0i(i 2n c 1 K( ( /)2,C2,t2|0i,ci,ti) = l (3.32) 

J d 2n <p 2 d 2n c 2 K(0 2 ,C 2 ,t2|01,Cl,tl) = 1 (3.33) 
and that it conserves the Grassmannian delta-function 5 2n (c): 



J d 2n <p x d 2n c x K(4>2,c 2 ,t 2 \4>i,ci,t 1 ) S^ l \ Cl ) = ^ 2n \c 2 ) 



(3.34) 



In fact, eq. (3.34) can be regarded as an expression of Liouville's theorem' 131 . Recall that 
the Liouville measure on phase space is given by 

dfx = cu n = ( 1 rf a A d<p b ) n (3.35) 
Consequently, if assumes its canonical, i.e., ^-independent form, we have 

= u n = n\ dcj) 1 A d(j) 2 A • ■ • A d<p 2n (3.36) 

Invoking the correspondence' 101 c a < > d(f) a between ghosts and differential forms on M.2n , 



* This is most easily seen by noting that M is of the form u> ab times a symmetric matrix 1151 and 
therefore lies in the lie-algebra of Sp(2n). 
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we see that 5 2n (c) corresponds to the Liouville measure 

6^ (c) = cW • • • c 2n <— > d/i (3.37) 

Eq. (3.34) expresses the fact that <i/i is invariant under the Hamiltonian flow. The 
infinitesimal version of eq. (3.34) is 

^(2")( c ) = (3.38) 
with the lie-derivative lh defined in (3.16). Eq. (3.38) follows from 

c^* (2n) (c)=0 (3.39) 

and d a h a = 0. It implies that not only zero-forms evolve according to the Liouville equa- 
tion (3.17) but also the coeflicent functions of the 2n-forms 

= Qi...2n{4>)c l c 2 ---c 2n EE g(<f>)S^ (c) (3.40) 

i.e., the equation dtQ = —Lg can be immediately derived from (3.38). In the follow- 
ing we shall use the 2n-form (3.40) in order to represent conventional scalar densities 
on phase-space. From a mathematical point of view the factor 5^ 2n \c) provides the vol- 
ume form on M.in , whereas from a physical point of view 5(c) is the vacuum of the 
fermionic Fock spaced In the Schroedinger-like picture mentioned before' 161 the condition 
c a \vac) = translates into c a < c\vac >= , i.e., the "position representation" of \vac) is 
< c\vac >= 5( 2n \c). 

To summarize this section we can say that in the classical path-integral the ghosts play a 
double role: a dynamical one, in the Heisenberg-like"'' picture of classical mechanics, because 
their equation of motion is the Jacobi equation and a geometrical one , in the Schroedinger- 
like picture of CM, because the time-independent c a span the cotangent space T£M.2n- 
This double role will be heavily exploited in our discussion of the Lyapunov exponents. The 
dynamical role of the ghosts had already been partly exploited in the second and the last 
of ref.[10]. There it was shown that for any hamiltonian H(<p) the action C is invariant 
under a set of universal (graded) symmetries which form an ISp(2) algebra. Part of this 



f See the appendix and ref.[16] for more details about this. 

\ We called it Heisenberg-like because the variables <p, c depend on t in this " picture" . 
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algebra consists of a BRS-like operator Q = ic a X a and an anti-BRS operator Q = 
icuj ab X b , respectively. £ was also shown to be invariant (up to surface terms) under a 
super symmetry generated by the charges 

Qh =c a (d a -d b H) 

3.41 

Q H =c a u ab (d b + d b H) 
which are nilpotent, Q H = Q H = and close on the superhamiltonian: 

in=[Q H ,Q H ] (3.42) 

In the second of ref.[10] it was shown that the phase of a dynamical hamiltonian system with 
this supersymmetry unbroken was the same as the ergodic phase of the system. Moreover in 
the last of ref.[10] the supersymmetry and its relation to ergodicity were applied in a study 
of the Toda criterion, which was the first criterion put forward in 1974 to detect transition 
from ordered to stochastic motion. The next step , in the physics-history^ of dynamical 
systems, was taken in 1974-75 and it consisted in using more refined tools to study these 
transitions. These tools were the KS-entropy and Lyapunov exponents to which we turn 
now. 



4. OBSERVABLES FOR HIGHER 
DIMENSIONAL LYAPUNOV EXPONENTS 

As we have already stressed the path-integral formulation' 101 of classical hamiltonian 
dynamics naturally involves the Jacobi fields c a (t). Therefore it seems plausible that it 
should be possible to relate quantities like the generalized Lyapunov exponents to certain 
observables in this theory. 

For any observable O = 0(<j), A, c, c) we define the "vacuum expectation value" as 

oo 

(0) = I V<f>VXDdDcO{<f>,\,c,c) 5 {2n \c(-oo)) exp i J dtZ (4.1) 

— oo 

where an integration over a (±oo) and c a (±oo) is understood^. Note that here we 
are dealing with a trace-formalism, rather than a bra-ket formalism. Therefore, contrary 

§ We say "physics-history" because in the mathematics literature (Kolmogoroff-Sinai ) these tools had 
been developed before. 

Here t=+oo (t=— oo) is a symbolic notation for some finite time which is larger (smaller) than any 
time argument of the fields in O. 
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to quantum mechanics, the "state", a generalized density g, appears only once under the 
path-integral' 101 . Expectation values of the type (4.1) can be reduced to strings of prop- 
agation kernels K(0j, cj, tj\<f>j-i, Cj-i, tj-i) connecting field monomials with different time 
arguments. An identity which we will often use is 

V<p VX Vc Vc A (0(t 2 ), X(t 2 ), c(t 2 ),c(t 2 )) ■ B(0(ti), A(ti), c(fi), c(h)) 



exp{i J dtC] e(0(ti),c(*i)) = 



ti 



(4.2) 



= I dfa dci A((f) 2 ,-i-^-,c 2 ,^-) 
J d(p 2 dc 2 



d d 

■ K(0 2 , c 2 , t 2 \(f)i, ci, ti)B(0i, -i-^— , ci, -z—)g(<i>i, ci) 

where A and B are arbitrary functions and where the functional integration is subject 
to the boundary conditions (f)(t 2 ) = <f) 2 and c(t 2 ) = c 2 fixed, while 4>{t\) and c(ti) are 
integrated over. For the observables considered in the following, A and B will be free from 
ordering ambiguities. Eq.(4.2) can be proven by discretizing the path-integral in the usual 
way. 

To start with we consider the family of observables 

O a b {T- O ) = c a (T)c b (0)S^ n Hm - O ) (4.3) 

for fixed T > and <pQ e M. 2n - The delta-function picks the trajectory <f)(t) which passes 
through the prescribed point <fio at time t=0. The creation operator q,(0) creates at t=0 a 
one-ghost state from the vacuum 5^ 2n \c) and c a (T) destroys it at some later time t=T. 
By using eq.(4.2) we obtain for the expectation value of O b : 



(O a b (T,<p )} = J #(oo)dc(oo)K(0(oo),c(oo),oo|0(T),c(T),T) 

• J d<j>(T) dc(T) c a (T) J d0(O)dc(O) K(0(T),c(T),T|0(O),c(O),O) 



J d(j){-oo)dc{-oo) K(0(O),c(O),O|0(-oo),c(-cx)),-cx)) 



• 5(0(0) - O ) (44) 

d 

• <5 (2n) (c(-oo)) 

The third factor of K on the RHS of (4.4) propagates the vacuum 5^ 2n \c{— oo)) 
from t = — oo to t — 0. Because of Liouville's theorem, eq.(3.34), the result at t — is 
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5( 2n )(c(0)). Then, between t=0 and t=T, the second factor of K propagates the ghost 
excitation created by ^(0) acting on the delta-function. For t > T we are left with the 
vacuum again. Using (3.32) we see that, again as a consequence of Liouville's theorem, the 
first K in (4.4) is ineffective. In this way (4.4) boils down to 

(0 & a (T,0 o )> = J d<t>{T) dc(T) c a (T) J de(0)K(<KT),c(T),T\<h,c(0),0)- 

° 5 (2n) (c(0)) (4.5) 



dc b (0) 

= J dc(T)dc(0) c a (T)5(c(T) - S(T; 0o)c(O))^<^)(c(O)) 

where also (3.31) has been used . Taking advantage of (3.39) the final result reads 

(O a b (T^o)) = S^T;<f )0 )9(T) (4.6) 

The step-function 6(T) has been included because (O) = for negative values of T, in 
this case in fact the destruction operator acts on the vacuum before the creation operator 
and therefore the probability is zero. We conclude that the two-point function of the ghosts 
is given by the Jacobi matrix S%. 

Later on we shall see that the higher dimensional Lyapunov exponents are related to 
observables of the form (0 < / < 2n) 

O f (T; O ) = c ai (T) • ■■c a f(T)c af (0) ■ ■■ c ai (0)5(0(0) - O ) (4.7) 

Their expectation values 

r f (T;<f> ) = (O f (T;4> )) ,T>0 (4.8) 

can be evaluated with the same method as above. The result is (up to an irrelevant factor) 

T f (T; O ) = S [a a \(T; O )^ 2 2 (T; O ) • • • S#(T; O ) (4.9) 

where the square brackets denote the complete antisymmetrization; for f=2, for example, 
(4.9) reads as T2 = (TrS) 2 — Tr(S 2 ). The interpretation of the observables Of is as 
follows. At time t=0 the operator c a/ (0) • • -c ai (0) creates a state with f ghosts (or better 
with f Jacobi-fields) from the vacuum. In the geometric interpretation of the theory this state 
corresponds to a f-form d<p ai A d<p a2 A • • • d<p af ,i.e., to a f-dimensional volume. Hence Tf 
contains information about the rate of growth of f-dimensional volume elements in tangent 
space. 
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The expectation values Tf are related to the Lyapunov exponents X\ (4>o) as follows. 
Consider first 

r 1 (t;<fo)=TrS(t;<l>o) (4.10) 

and assume that the initial point 0o gives rise to a periodic trajectory with period r,i.e., 
4> a (t) = 4> a (t + r). Consequently the matrix M of eq.(3.28) is periodic too, and Floquet 

[171 

theory tells us that S can be written as 

S(t) = P(t)exp(Rt) (4.11) 

where P is a periodic matrix, P(t) = P(t + r), and R is a constant one. Because 
5(0) = P(0) = 1 we have S(r) = exp(Rr). Let us then diagonalize R; its eigenval- 
ues p a , a = 1 • • • 2n, are the characteristic exponents of M{t) and exp(p a t) the correspond- 
ing Floquet multipliers. In the eigenvector basis of R we have (suppressing the argument 
0o) 

S£(t) = P b a (t)exp(p a t) (4.12) 

We assume the ordering 

Repi > Rep 2 >■■■> Rep 2n (4.13) 
If Repi is strictly larger than Rep 2 , the large-t behaviour of Ti is 

ri(t;0o) ~P\exp( Pl t) (4.14) 

Hence 

01 (0o) = limsup -lnFi(t; <po) = pi (4.15) 

exists and coincides with the Lyapunov exponent A^(0o)- 

Because S(t) is a real symplectic matrix, its eigenvalues appear as 4-tuples: if p is 
a complex eigenvalue then' 141 ^ and are eigenvalues too (not necessarily different 
from p). Since S(r) = exp(Rr) this means that if p is a characteristic exponent, then 
also p*,—p and — p* are characteristic exponents. Let us write p a = l a + iuj a with 
l a and uj a real. As for the relative magnitude of the numbers Rep a = l a two cases 
have to be distinguished. If some p a has a nonvanishing imaginary part, uj a ^ 0, then 
p* is different from p a and consequently there exist two exponents with the same real part 
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(they give rise to equal contributions to the Lyapunov exponents). On the other hand, if 
uj a = , the exponent p a is real and generically there will be no other exponent with the 
same real part. In this case Ti, say, is dominated by a single real multiplier, i.e., pi. If 
however, p\ = l\ + ioJ\ and P2 = h — iu)\ are a complex conjugate pair, eq.(4.14) is 
replaced by 

ri(f; O ) ~[R\(t) e iuJlt + R 2 2 e~ luJlt ]e ht (4.16) 

From (4.9) with (4.12) we obtain for the higher correlation functions 

T f (t; O ) = P [a a\{t) ■ ■ ■ P£l\t) exp[(p ai +Pa2 + ... + Paf )t] (4.17) 

Due to the antisymmetrization, the indices a,j of the p's in the exponential must all be 
different. Because of the ordering (4.13) this implies that for t — > oo 

T f (t; O ) ~ p(t) exp[( Pl + P2 + --- + p f )t] (4.18) 

for some r-periodic function p(t). As it stands (4.18) is correct only if the real part 
of the last eigenvalue, If = Repj, is strictly larger than the real part of the following 
eigenvalue, lf+i- If, for some reason If = lf+\, the asymptotic formula consists of two 
terms, 

T/(t;0o) exp[(pi + p 2 H + Pf)t] + 

(4.19) 

p(t) exp [(pi + P2 + --- + Pf-i + Pf+i)t] 
or even more terms if pf is degenerate with p/+2; P/+3 > ' ' ' ; e ^ c - 111 an Y case 




is the sum of the f largest real parts of the eigenvalues p a . 

In the introduction we mentioned already that the leading Lyapunov exponent governing 
the evolution of f-forms, X[ (<fio), is related to the higher Lyapunov exponents for one forms, 
A^((/>o), according to 

/ 

x[ f) (M = J2 X ^^) (4.21) 

i=i 

This is exactly the relation we found in eq. (4.20). The correlation function is the expec- 
tation value of the operator Of which creates and destroys a f-dimensional "parallelotope", 

* e.g., because p/ and form a complex conjugate pair. 
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and hence gj describes the rate of exponential growth of f-dimensional volume elements 
in tangent space: A^(</>o) = 9f(<fto)- Once gj is known for all values of f, the system of 
equations (4.20) can be solved for k((/)o) = X^i^o) in order to obtain the higher Lyapunov 
exponents for one forms, i.e., A^(0o)- 

So far we have shown that the correlation functions (0(t; (fio)) encapsulate the informa- 
tion about all the Lyapunov exponents related to a fixed trajectory, namely the one starting 
at <pQ a time t=0. If we restrict <fio to a region in phase-space of connected stochasticity 
(excluding regions of regular motion) the Lyapunov exponents are independent of 0o and we 
may equally well extract them from the observable (4.7) with the delta-function fixing the 
initial point omitted, but under a path-integral which is over a restricted class of trajectories 
only. (As was shown by Oseledec 1131 it is not really necessary to insist on closed trajectories.) 
As we shall see in the next section, expectation values of this type are closely related to the 
partititon function of the superhamiltonian 7i. 

5. PARTITION FUNCTIONS OF THE SUPERHAMILTONIAN 

Obviously the superhamiltonian H of eqs. (3.15), (3.16) does not mix forms of different 
degree (ghost number). Therefore it makes sense to consider H p , the restriction of H to the 
space of homogeneous p-forms^, which is spanned by the generalized densities of the type 
(3.8). Let us evaluate the partition function 

Z P (T) = Tr[exp(-mt)} = £ < Xp*\exp(-iH p t)\ Xp * > (5.1) 

a 

where is a basis of p-forms and {(x p |} the dual basis of p- vectors (antisymmetric 

contravariant tensors of degree p). In component notation the completness relation reads 

Ex^w-a x^'y bi '" bp = 6(2n) ^ - ■ ( 5 - 2 ) 

a 

From (3.7) and (3.31) we obtain for the matrix element of exp(—iH p t) 
(0, c\ exp(-iH p t ) \ X p) =< ^ d \t; 0), S^ 1 ^; 0)c|Xp >= 

= ^(^(tl <P)) ai ... ap S~\t; </>)% ■ ■ ■ S-\t; ...^ (5 ' 3) 

The partition function (5.1) is obtained by stripping off the ghosts from (5.3) , contracting 
with the dual basis and integrating over <ft . Exploiting (5.2) one finds (up to an unimportant 

| For an analogous construction in supersymmetric quantum mechanics see ref[18]. 
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constant) 

Z P (T) = [ dcf>5(<f>-$-\T;cf>))S-\T; ( f>)t 1 ---S-\T;ci>) a a p J 

~ J d<t>5{® d {T-4>)-<t>) S(T;4>)t---S(T-4>) a a T n Z P J 

In the second line we used the identity [A. 55] (which is derived in the appendix) and the 
fact that the Jacobi matrix is unimodular. As it was to be expected, Z P (T) receives 
contributions only from closed trajectories of period T. 

Next we show that the partition functions in the p-form and the (2n-p)-form sector 
coincide: 

Z p (T) = Z 2n - P {T) (5.5) 

The reason is that there exists a duality operation * which maps p-forms on (2n-p)-forms 
and which commutes with H. On an arbitrary p-form \ it ac ts as 

(*x)a P+1 ...a 2n ~ e ai ... a2n u a ^ ■ ■ • oM> Xbi ... 6p (5.6) 

This kind of * operator is analogous to the Hodge operator on Riemannian manifolds. The 
★-operation commutes with 7i = —ilh because the Lie-derivative along the hamiltonian 
vector field of both the e-tensor and of the uo ah vanishes. (The equation lh^ai—a p — is 
the component form of eq. (3.38).) Hence the spectra of 7i p and H,2 n -p coincide, which is 
analogous to the well-known Poincare duality for the Laplacian. From (5.4) with (5.5) and 
(4.9) it follows that 

Z f (T) = fdcf> 5($ cl (T; 0) - 0) S(T; <f>) [a a \ ■ ■ ■ S(T; 4>)V) 

(5.7) 

# 5($ d (T;0)-0) r7(T;0) 



Thus we find that the ratio z f ) T ) can be interpreted as the average of Tj(T) over all closed 
trajectories of (not necessarily primitive) period T. Because Tf={Of) 
with Of given in (4.7) , we easily could write down a path- integral representation for 
Zf by combining (5.7) with (4.1). Note also that 

Z (T) = J ^5($ d (T;0)-0) (5.8) 

"counts" the number of initial points which lead to a closed trajectory of period T. Of 
course these points are infinite in number. What we should do in (5.8) , and in all the 
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other Zf(T), is to factor out the equivalence relation which relates two points which are on 
the same trajectory. This work is in progress' 191 . After having factored out the equivalence 
relation, the new Zq(T) would count the number of closed orbits of period T and that implies 
that from this new Zq(T) we should be able to get the topological entropy of the system' 201 . 
Anyhow the reader should not be worried about the "infinity" produced by Z P (T) because 
the physical quantitities which we will consider are always given by ratios of Z P {T) and in 
these ratios the infinity cancels out. 

Until now we were purposely vague about the domain of the <fi integration in eq. (5.7). 
If the trace in (5.1) is over p-forms defined on the full phase-space M.2n then clearly the 
integration in (5.7) is over all of M.2n- However, in general one would like to discuss 
the chaoticity properties of a system in different regions of phase-space separately, and, 
most importantly one would like to work at fixed energy E. The (2n-l)-dimensional energy 
hypersurface M.2n-i{E) is the subspace of M.2n on which H(4>) = E. If the initial 
condition of the classical path-integral are fixed such that the initial point 0(0) lies on 
M.2n-i {E), then the dynamics is such that (p(t > 0) is still on M.2n- Correspondingly, 
if a zero-form g(<p, t = 0) has support on M.2n-l{E) only, this property is conserved 
under the time evolution. This is not sufficient, however. We also have to make sure that 
p-forms on A42n-l(E) evolve into p-forms on M.2n-\ (E). A p-form on M.2n-l(E) is a 
tensor which has no components in the direction perpendicular to the energy hypersurface. 
Loosely speaking, the exterior algebra on M.2n-i{E) is obtained by putting to zero the 
component of d(f) a normal to the energy hypersurface: d a H(4>)d(f) a = 0. Therefore, defining 
N(t) = d a H((J)(t))c a (t), we have to constrain the path-integration to the subspace with 
N(t) = 0. In the second of ref.[10] we have shown that the charge N is conserved under 
the time evolution: [iV, H] = 0. In fact, N is the difference between the supersymmetry 
generator and the BRS-generator: N = Qh~Q- This implies that, imposing N(t = 0) = 0, 
guarantees that N(t)=0 at any later time. As a consequence, if g((f),c,t = 0) is a tensor 
on M.2n~i{E), (i.e., if it does not contain any factor of N), also g((f),c,t) at t > is 
a tensor on M.2 n -i(E). Thus we can consistently truncate the classical path-integral to 
the energy hypersurface. In particular we may define partition functions Z p (T; E) as in 
eq.(5.1) but with x P a complete set of p-forms on M.2n-i{E). The next step is to develop 
a full symplectic and coordinate free formalism. This implies that we will have to decrease 
the dimension of M.2n-i{E) by one unit to go to an even-dimensional subspace of M.2n ■ 

[21] 

This is done via the so-called Batalin-Fradkin-Vilkovisky method (BFV-formalism) which 
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implies the introduction of further auxiliary fields and further ghosts. The details of this 

[191 

will be presented elsewhere 

Let us now go back to Zf(T). From the asymptotic T — > oo behaviour of Zf(T) we 
define the generalized Lyapunov exponents Ay, 1 < / < 2n , according to 



Z f (T) 



exp[(A 1 + A 2 + --- + A f )T] (5.9) 



Zo(T) 

This is the deterministic analogue of eq.(A.57) for stochastic systems which we will find 
in appendix A. Formally the exponent Ai, say, is defined as in eq.(2.11) with the only 
difference that the ensemble average (•) is not taken with respect to the stochastic measure 
but with the deterministic one. Because we formulated also the deterministic systems in 
a path-integral language, this correspondence becomes particularly transparent. Note that 
generically (for a deterministic system) there is no simple relation between the generalized 
exponents Aj and the ordinary ones, A^. In order to obtain the former, one averages the 
monodromy matrix S(T) for many paths of length T and sends T to infinity afterwards, 
whereas the latter one is obtained from the large-T behaviour of S(T) on a single trajectory. 

We now briefly comment on the eigenvalue problem of the superhamiltonian 7i and the 
time evolution operator exp(—i7it) . It is at this point that we encounter the most important 
differences between stochastic systems and classical hamiltonian systems. In the former case 
(see appendix A) the superhamiltonian is a second order Schroedinger operator, in the latter 
it is the first order Lie-derivative operator H = —ilh- Let c ) be an eigenfunction of 
H in the p-form sector so that 

exp{— iHt)x{4>, c) = exp{— i£t)x(4>, c) (5.10) 

for some constant £. Using eq. (5.3) we immediately see that the components Xai-a p {<t>) sat- 
isfy the eq. 

Xai -a p {<$>cl(t; <P))S(t; 4>)H ■ ■ ■ S(t; 4>)H = exp(+iSt) Xbl -bM ( 5 - U ) 

On the LHS of this equation we recognize the usual tensorial transformation law under the 
Hamiltonian flow. It affects the eigenfunctions x only via the overall factor exp(iSt). Let 
us look at p=l in more detail, where the relation is 

Xa(*d(t; <P))S(t; 4>) a b = exp(i£)xM (5.12) 

and it has to hold for all <fi and all t if x is an eigenfunction. Let us assume we pick a point 
00 on ■M.2n which is the initial point of a closed trajectory of duration r: $ c ;(r; 0o) — 00- 
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Let S(t; 0o) be the the Jacobi matrix evaluated at t = r*. If we evaluate (5.12) for the 
special values <p = <po and t = r , we find 

[5(r; 0o)b - e^(^r)5 & a ]xa( ( />o) = (5.13) 

i.e., that exp{i£r) is an eigenvalue of the monodromy matrix provided that Xa(4>o) ^ 0. 
This shows that the possible values of £ are closely related to the periods and the Floquet 
multipliers of closed orbits. 

As a final remark we mention that in a previous paper (the fourth of ref.[10]) we had 
shown that the alternating sum 

2n 

Z(T) = STr[exp{-iHt)] = J^(-1)%(T) (5.14) 

is a topological invariant (the Euler number of M.2n) and has no dynamical significance 
therefore. It was proven in fact that the alternating sum (5.14) is invariant under deforma- 
tion of the hamiltonian vector fields. On the other hand, we saw that the individual ZpS 
are not invariant under such deformations and therefore can be used to characterize certain 
properties of the dynamics. 

6. CONCLUSIONS 

To summarize we can say that, looking back at eqs. (4.7), (4.8), (4.20)for the ordinary 
Lyapunov exponents and at (5.9) and (5.1)for the generalized ones, both exponents admit a 
very simple and natural representation in terms of classical path-integrals. The representation 
of the Lyapunov exponents as expectation values of some observables allows, for example, 
then for a perturbative calculations of them in the same manner as it is usually done in 
field-theory. Viceversa the representation via eq. (5.9) allows the use of spectral methods 
to calculate then. The supersymmetry of the path-integral is crucial in this context: it re- 
lates the "bosonic" dynamics of the trajectories 4>(t) to the evolution of its "fermionic" 
superpartner, the Jacobi-field c(t) and because of this relation it was natural to expect 
that information on the dynamics of the Jacobi-fields could be extracted from some objects 
containing only the dynamics of the standard-phase-space variables as Zf is. This super- 
symmetry will also produce' 161 semplifications in the perturbative calculations in the same 
way as it does in standard field theory. 

[22] 

* This is what is called in the literature "monodromy matrix" of this loop. 
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A lot of work remains to be done in order to extract the dependence on the energy of the 
ordinary (and generalized) Lyapunov exponents and of the various entropy-like quantities 

[19] 

that are associated to them. This work is in progress together with an understanding of 

[23] 

the relation of our formalism with the thermodynamic formalism of Ruelle on which we 
briefly comment in appendix B. 



APPENDIX A 

In this appendix we use the tools of supersymmetric quantum mechanics in order to dis- 
cuss the generalized Lyapunov exponents for stochastic systems. The reader should compare 
the various steps of the derivation with their classical counterparts described in the main 
body of the paper. 

Let us consider a stochastic process on an N-dimensional, metrically flat configuration 
space with (local) coordinates X{ i = 1,---,N. The dynamics of the random variables 
Xi(t) is given by the Langevin equation 

ii(t) = -diU(x(t)) +rH(t) (A.l) 

where U(x) is a smooth potential and rji(t) is a white noise : 

(»*(*)) = o, <»*(t)i&(0> = M(*-0 ( A - 2 ) 

ryi 

It is well known that the stochastic correlations 

{x v ■ ■ ■ x v ) derived from (A.l) can be obtained from a supersymmetric generating func- 
tional of the form: 

Zsusy = J VxV^jViJj exp(— J dtLsusy + source terms) (A.3) 

where 

1 1 

Lsusy = -xf + -(diU) 2 + fa [dtdij + didjUix)]^ (A.4) 

The Xi are commuting variables while ipi and ipi are Grassmannian ones. The supersym- 
metry transformations under which (A.4) is invariant are given by: 

dxi = - eipi + etpi 

Sfo =e(-Xi + diU) (A.5) 
54>i = e(xi + dill) 

In a Schroedinger picture formulation of the supersymmetric quantum 
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mechanics' 16 ' 181 denned by (A. 4) the states |<&) are described by wave functions 

$(xi,i/>i)=<Xi,il)i\$> (A.6) 

depending on the variables xi,ipi. The operators Xi and vf)i acton $(x,ift) by multipli- 
cation and their conjugate momenta by differentiation: 

d d 

Eq.(A.4) gives rise to the (Weyl ordered) superhamiltonian: 

Hsusy = Hb + Hp (A. 8) 



where the "bosonic" part is: 



H B = -\d 2 + ^diUdiU, d 2 = (A.9) 



and the "fermionic" part is: 



= \d 2 U{x) - didjUix)^ — 



In terms of the supercharges 



Q =[di + diU{x)} 



di/ji 



we have 



2 

[16] 



(A.10) 



Q = [-d l + d i U(x)]tp i 

d (A.ll) 



Hsusy =\[Q,Q] (A.12) 



In this operatorial formalism a generic wave function $(x, ip) possesses an expansion 
of the form 

$ (^) = EAW fcl -^ (A.13) 



which is reminiscent of the expansion of an inhomogeneous differential form in a basis 
dx kl A • • • A dx kq . We say that $ has ghost number "p" if on the RHS of (A.13) only the 
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term with q=p is different from zero (homogeneous form of degree p): 

V) = -M.-k, W kl • • • (A.14) 
pi 

The vacuum of the fermionic Fock space, defined by vpi \ vac) = 0, is represented by a wave 
function of ghost number n: 

< x,ip\vac>=$ vac (x)8(ip) (A. 15) 
where the Grassmannian delta-function 

Sty) =^^ 2 ...^ N (A.16) 



can be visualized as describing a completely filled "Dirac sea". Multiparticle states are 

< _ _d_ 



obtained from \vac) by acting on Sty) with the "creation operator" ip i = -Mr-. A state 



containing f "particles" has a wave function with ghost number N-f: 

<x ,W f >=**... kf (x)^.--^6M (A.17) 

The time evolution of the states ip, t) is governed by the Schroedinger equation 

1 d d 

H(x, - — qtMx, V>, t) = -d t <f>(x, V, t) (A.18) 

It has the formal solution (t > 0) 

ip,t) = J d N x d N ip K(x, ip, t\x , ip , t )<f>(x ,ip , to) (A.19) 

The evolution kernel K is a solution of the Schroedinger equation (A.18) with initial 
condition 

K(x, ^ toko, Vt), t ) = 8 N (x- x )8 N ty - iPo) (A.20) 
Its path integral representation involves the lagrangian (A. 4): 

K(x 2 ,^2,t2|xi,^i,ti) = J Vx(t) Vip(t) V$(t) exp{- J dtL susy ) (A.21) 

The boundary conditions are x{t\^) = £1,2, ^(^1,2) — ipi,2 an d ^(^1,2) is integrated over. 
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\{x 2 ^ 2 MWMM) = Jvx{t) exp{- J dt[^a$ + ^(diU) 2 ]y 



(A.22) 



(For details see ref. [16,12]) Writing the above kernel as 

h 

U 

■ K F (^2,t 2 \ipi,tr, [x]) 

with 

K F (^ 2 ,t 2 |^i,ti; [x]) = J Vi> V4> exp{- J dt^i[dtSij + didjU{x)]ipj} (A.23) 

we can explicitly evaluate the fermionic kernel Kf which is a functional of the bosonic path 
X{(t). Following the treatment of ref. [12], we first perform the (unconstrained) if) -integration 
in eq.(A.23) which leads to 

K F ftM2|^Mi; W) = SiidtSij + didjUix))^} (A.24) 

Obviously only the solutions of the equation 

fa = -didjU{x{t))^j (A.25) 

obeying the boundary conditions ^(^1,2) — i>l,2 contribute to Kp . This is a rather 
remarkable fact, because eq.(A.25) is precisely the Jacobi eq. pertaining to the Langevin 
equation (A.l) (i.e., if Xi(t) is a solution of (A.l), then the variation 8xi(t) = ifti(t) is a 
solution of (A.25)). The explicit solution of eq.(A.25) reads 

tl H (t) = S lJ (t;[x})^ j (t 1 ) (A.26) 

with the Jacobi matrix 

t 

S(t; [x]) = fexp{- J dt'M(x(t'))} (A.27) 

(T denotes the time ordering operator) and where 

Mij(x(t))=didjU(x(t)) (A.28) 

S is a solution of the matrix equation S = —MS with S(ti) = 1. We obtain for the kernel 
(A.24) 

K F (V>2 , t 2 \ipi,ti; [x]) = 

6[l>(t) - S(t; [x])1> tl ] det[d t 5ij + d l d j U(x)] (A ' 29) 

In order to give a well-defined meaning to the determinant, we have to specify a discretization 
scheme for the functional integral. Because the Hamiltonian in (A. 8) was chosen to be Weyl 
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ordered, we must use the mid point rule for the discretization.' 241 For this discretization, it 



is known that [25] 



12 

det[d t 5ij + didjU] = exp{^ J dtd 2 U(x(t))} (A.30) 

Hence the final result for the fermionic kernel is 

K F (ip 2 ,t 2 | ipi,t\\ [x]) = 

ti 

= 6(ij 2 - S(t 2 ; [x])ih)exp{^ J dtd 2 U(x(t))} (A ' 31) 

It is possible to check eq.(A.31) also without referring to path-integral manipulations. The 
path-integral on the RHS of eq. (A. 23) is the formal solution of the Schroedinger equation: 

[\dfu(x) - didjU{x)^~ + fljKH^M^Mn W) = o (A.32) 

with Ki?('0,£i|'i/;i,ii; [x\) = 5(ip — ipi). It can be checked that the RHS of (A. 31) does 
indeed solve this initial value problem (for t > 0). For future reference we note that Kp can 
also be written as 



K F (ip 2 ,t 2 | ipi,h; [x]) = 



(A.33) 



= 5{S(t 2 ;[x])- l ^ 2 -^ 1 )exp{~ J dtd 2 U(x(t))} 

ti 

because S = —MS implies that 

t 

det[S(t)] =exp{- J dt'd 2 U(x(t')} (A.34) 
ti 

Obviously the Hilbert space Ti of our theory is a direct sum of subspaces with a fixed ghost 
number p, < p < N. Instead of labelling the various sectors by their ghost number (or, 
equivalently, their degree as a differential form) we shall use f = N — p, < / < N, so 
that the subspace Tif of H. consists of wave functions of the form (A. 17). Consequently 

N 

«=©«/ (A-35) 
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The partition function decomposes correspondingly 

ZsusyiT) = Tr[exp(-TH susy )\ 

A , N (A.36) 

with 

Zf(T) = £ < $ /l exp(-TH f ) \<S>J) (A.37) 

where Hf denotes the restriction of H susy to TCf, and is a complete set of states 
with ghost number p=N-f. The Hamiltonians are particularly simple for f=0 and f=N. In 
these cases they are diagonal in the tensor indices: 



Ho = - -d 2 + -diUdiU - -d 2 U 

I I I (A.38) 

H N = - -d 2 + -diUd t U + -d 2 U 

The respective partition functions can be represented by purely bosonic path-integrals: 

T 

Zq,n( t ) = J Vx ( l ) ex P{~ J dt L o,n} (A.39) 

pbc 

where 



Lo = -x 2 + -d l Ud l U--d 2 U 

1 9 1 „ . 1 



(A.40) 

L N = ^-x 2 + ^-diUdiU + ^-d 2 U 



2 ' 2 2 

are the restriction of L susy to Ho and Hn, respectively. In (A.39) we used periodic 
boundary conditions (pbc): x(0) = x(T). As a preparation for the discussion of the Lya- 
punov exponents, we shall now give an alternative representation of Zo and Zf in terms 
of the complete supersymmetric lagrangian L susy of (A. 4). For f=0 we may write 

r 

Zo(T) = J Vx(t) J ViP(t)V4>(t) exp{- J dt L susy }5(iP(0)) (A.41) 

pbc 

Again, periodic boundary conditions are used for the Xj-integration, but if)(0) and ifi{T) are 
treated as independent integration variables. Note that these boundary conditions are dif- 
ferent from the ones usually used in supersymmetric quantum mechanics. The equivalence 
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of (A. 41) and (A. 39) is established by noting that: 

T 

Jvi/j exp{- J ^[dtSij +didjU]^j}5{^)) = 
o 

= J #(T) #(0) K F (^(T),T^(0),0; [x])5(^(0)) (A.42) 

T 

= ex v{\J dtd 2 U{x{t))} 
o 

where (A. 31) has been used in the last line. The term ~ d 2 U found in (A.42) , when 
combined with the first two terms on the RHS of eq.(A.4), yields exactly the lagrangian 
Lq of (A. 40). In the same way one can show that 

T 

Z N (T) = J Vx(t) J ViP(t)Vi>(t) 8(i;(T)) exp{- J dtL susy ) (A.43) 

pbc 

(Here and in the following we ignore multiplicative constants.) 

In deriving eq. (A. 31) we realized already that the dynamics of the Grassmannian 
variables ipi(t) is essentially deterministic. In fact, given a fixed trajectory Xi(t), the kernel 
Kf is non-zero only if the initial and final value of if) are related by the " classical" Jacobi 
matrix Sfe; [x]). Because all the information about Lyapunov exponents is contained in 
S, this suggests that it should be possible to extract them from certain correlation functions 
involving ghosts and antighosts as we did in the deterministic case. In fact, let us consider 

F f (T) = <^ 1 (T)^ 2 (T)---^ / (T)^ / (0)-- -^(0)^(0)) (A.44) 

with the expectation value (O) of observables O defined by 

T 

(O) = ZoiT)- 1 J Vx J V^ViP O 8(iP(0)) exp{- J dtL susy ) (A.45) 

pbc 

where ip(T) and ^(0) are independent integration variables. As we discussed in section 
4, the operator ijj] if ■ ■ ■ i/j^ creates at time t=0 a f-volume from the vacuum which is 
propagated by the supersymmetric dynamics until it is destroyed by ip^ ■ ■ ■ vpk f at time 
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t=T. This becomes obvious if we use (A. 4) in (A. 45), in order to write 

T 

T f (T) = Z 1 J Vxexp{-J dt[h 2 l+ \{diUf\}Qf{TM) (A.46) 

pbc 

with 

T 

Gf(T; [x]) = J Vi) V^P exp{- J dt ^ + didjU]^} 

o 

• rf, h (T) • • • ^ fc/ (T)4, kf (0) • • • ^ (0) 6M0)) 
= J d N tp(T) d N m ipkAT) ■ ■ -ipk f (T) K F (^(T),T\m,0;[x}) 

mo)) 



dip kf (0) dip kl (0) 
Using (A. 31) it is easy to abtain: 



T 



Q f (T, [x]) = exp{\ J dt d 2 U(x)} ■ S^S k k l ■ ■ ■ S k k f (A.48) 
o 

with Sj = Sij(T, [x]) (we do not distinguish between upper and lower indices here). Inserting 
(A.48) into (A.46), and making use of (A. 40) , we obtain the final result 

r/(T) = ((4X 2 --< /] )} (A.49) 
with the average ((• • •)) performed by means of the lagrangian Lq, i.e., 

T 

({0[x})) = Zq(T) _1 J Vx 0[x] exp{- J dtL } (A.50) 

pbc 

Before relating (A. 49) to the Lyapunov exponents, let us briefly discuss another represen- 
tation of Ff(T). We shall show that Tf is the ratio of the partition functions of H and 
Hq, respectively: 

Z f (T)_ Tr[ex P (-TH f )} 
/( )_ Z (0) ~ Tr[exp(-TH )] ^ 

In order to evaluate Zf as written down in eq. (A. 37) , we first note that the components 
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of the vector exp(—THf)§J can be represented by the bosonic functional integral 

T 

(exp(-TH f )$])(x,ij) = J Vx exp{- J dtL N }$ a f (x(Q),S-\T,[x\)ip) (A.52) 

o 

with x(T) = x and where Vx includes an integration over x(0). To arrive at eq. (A.52) , 
we used (A. 20) with (A. 22) and (A. 33), as well as the definition of Ln in (A. 40). 
Furthermore, the inner product of two p- forms of the type (A. 14) is given by 

< ^ (1) l^ (2) >=jj d N x (A.53) 

If we use (A.52) with (A.53) in (A. 37), and exploit the completness relation of the we 
arrive at: 

T 

Zf(T) = jvx exp{- J • " • (S-'tli ( A - 54 ) 

pbc 

The occurence of N-f factors of the Jacobi matrix is due to the fact that Tij consists of 
moniomials with N-f factors of ip, i.e., differential forms of degree N-f. However, if we use 
the identity 

4" 1 • • • $ = 4 det(S)(S~^ ■ ■ ■ (S-'ti (A-55) 

with det(S) given by eq. (A. 34), and Cjy some constants, we see that Zf(T) coincides with 
Z Q (T) ■ T f (T) as given in eq. (A.49) . This completes the proof of (A.51). 

Benzi et al. [81 and Graham' 91 have studied eq. (A.49) for the special case of a one- 
dimensional configuration space, i.e., N=l. They determined a generalized Lyapunov expo 
nent Ai from the behaviour of T\(T) for T — > oo : Ti(T) ~ exp{K\T). If we represent 
Ti as the ratio of two partition functions, like in eq. (A.51), we see that Ai can be expressed 
in terms of the lowest eigenvalues E™ m and E™ m of Hq and Hi , respectively: 

Ai = E^ in - Ef n (A.56) 

If supersymmetry is unbroken, the vacuum is non-degenerate so that 
either E™ m < E™ in or Eq 1171 > E™ m . Only in the second case, when the vacuum is 
in the one-form sector, a positive Lyapunov exponent is obtained. 
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In the work of Benzi et al. and Graham only one dimensional systems were con- 
sidered, so that only the exponent Ai related to the evolution of one-dimensional volume 
elements could be defined. In the present paper we suggest the following interpretation of 
the correlation function Tf(T) of (A. 49). The antisymmetrized product of S- matrices de- 
scribes the evolution of f-dimensional volume elements. Therefore the large-T behaviour of 
Tf(T) can be used to define higher dimensional Lyapunov-like exponents: Ay, 1 < / < N. 
In view of eq.(2.10) and the discussion in section 4, it is natural to define 

F f (T) ~exp[(A 1 + A 2 + ---Af)T] (A.57) 

for T — > oo. This is the higher dimensional generalization of eq. (2.11) discussed in the 
introduction. By virtue of eq. (A. 51) the sum of the first f A's is given by the lowest 
eigenvalue E ! J J%n in the f-form sector: 

Ai + A 2 + • • • + A f = E^ in - Ef n (A.58) 

The above relation was conjectured in ref.[9] but no proof was provided, while here, using 
the machinery developed above, we have given a proof of it. For the deterministic case we 
cannot have a relation like (A.58) because the spectrum of the relative Hamiltonian TLf is 
not bounded below 
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APPENDIX B 



In this appendix we briefly comment on the relation between our path-integral formal- 

[23] 

ism and the thermodynamic one of Ruelle . The basic object of ref.[23] is a Fredholm 
determinant that is by now called the Ruelle Zeta-function ((z) is defined as 

oo n 

C(z) = det[l - zK {0) ] = exp [- £ Z - Tr(Kfa)] (B.l) 

n—l 

where z is a complex variable and Ko (for the classical Hamiltonian systems) is the kernel 
represented in eq. (3.9) with the ghosts c a put to zero and the interval of time taken to be 
a finite one which we will call A. So, roughly speaking, K(A) are the matrix elements 
of exp [— iHoA] which is the operator of evolution for zero-forms. The above series can 
be defined for a wide class of maps and it has a finite radius of convergence. The operator 
exp [— iHoA] is also called the Koopman operator and its inverse the Perron-Frobenius 
operator 1261 Using the formalism of Ruelle, it is easy to prove that the series (B.l) is also 
equal to the following one 

oo n 

C( z ) = exp[-J2-an] (B.2) 
z — ' n 

n=l 

where, (for the Hamiltonian evolution of zero-forms), 

a n = number of periodic trajectories of period nA (B.3) 

Note that, using eqn. (5.8), we can say that 

a n = Z {nA) (B.4) 

Here Zq has to be evaluated with all the cares we indicated under equations (5.8), basically 
following the lines of the last of refs.[23]. Ruelle proved the relation (B.2) for very general 
maps and to do that he had to use a highly sophisticated mathematical machinery. At a 
formal level, for classical hamiltonian systems, the proof is straightforward if we use our 
path-integral representation' 101 . In fact, following the steps (3.4) through (3.9) we see that 
the trace of the Koopman operator boils down to an integral over a Dirac delta function 

Tr[exp(-iH (nA))] = Z (nA) 

and this proves the theorem above. We do not pretend any mathematical rigorosity in the 
proof we have given above, but we have provided it to give to the reader some intuition 
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about the Ruelle Zeta-function. Another point to be cautious in the above proof is the 
problem of the continuum limit. The Ruelle zeta function was always built for discrete 
maps and not for the continuous time evolution given by our path-integral. So in the proof 
above, Zq(uA) should be considered the finite-lattice approximation to our path-integral, 
and actually that is how path-integrals are defined. Taking the limit of A going to zero has 

[27] 

to be done with the same care as explained in the book of Feynman-Hibbs for standard 
path-integral. If we want to avoid that, then we could choose A as long as the whole interval 
of time we are interested. Doing that we can then use the Ruelle Zeta-function to derive our 
Zq(A), in fact it is easy to see that 

Z (A) = Tr[exp (-i(W A))] =-f^*M (B.5) 

So one see that, using the appropriate zeta-function, we can extract one of the partition 
functions we have used in the paper. 

The same construction can also be applied to the kernel of evolution of higher forms and 

[23] 

this was also already envisioned by Ruelle in 1976 . We can define a generalized Ruelle 
zeta-function as 

oo n 

C(p)(*) = det[l - zK {p) ] = exp [- £ ?-Tr(Kfa] (B.6) 

n=l 

where we called K/ p \ the operator exp — %H V A appearing in (5.1) which makes the evolution 
of p-forms. Like for the zero-form case, for these generalized zeta functions (z) there exists 
an alternative representation, namely 

Tl 

C {p) {z) = exp[-Y, Z -o% ] ] (B.7) 

(p) 

where the coefficients are given by 

a n ] = E <P P er) [ al • • • S~\nA- ^rfaj] (B.8) 

Here S% is the Jacobi matrix we introduced in eq. (2.3) and 4> per denotes all points from 
which a periodic orbit of period nA originates. From eq. (B.8) and (5.4) we get immediately 
that 

ajf } = Z p (nA) (B.9) 

The proof of the above theorem is basically what is contained in formulas (5.1) through (5.4). 
One could have also used directly the K of eq. (3.9) and have all the higher forms included 
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in a single formula: 

oo 

C(z) = det[l - zK] = exp[- ^ Tr(K n )] (B.10) 
The trace TrK n can be taken either by choosing periodic boundary conditions or antiperi- 

[23] 

odic ones for the ghosts. In the second case, which corresponds to Ruelle's choice , we 
will get an expression which is an alternating sum (according to the order of the forms) 
and which will not depend at all on the Hamiltonian. This has been proved in the third of 

[28] 

refs.[10] along the lines of modern topological field theory . It is remarkablethat the same 
trick of getting something which does not depend on the Hamiltonian, by working with an 
alternating sum, was already employed by Ruelle in 1976. In our paper on the contrary, we 
have exploited also the individual terms of this alternating sum, and seen that they are re- 
lated to higher order Lyapunov exponents. Finally we mention also the work of Christiansen 

[29] 

et al. where forms were used in order to get information on the dynamics and not just on 
topological features of the manifold. In this work' 291 forms have a different meaning than in 
our formalism, however. 
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